On the Glauber model in a quantum representation 
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Abstract 

The Glauber model is reconsidered based on a quantum formulation of the Master equation. 
Unlike the conventional approach the temperature and the Ising energy are included from the 
beginning by introducing a Heisenberg-like picture of the second quantized operators. This method 
enables us to get an exact expression for the transition rate of a single flip-process Wi{ai) which 
is in accordance with the principle of detailed balance. The transition rate differs significantly 
from the conventional one due to Glauber in the low temperature regime. Here the behavior is 
controlled by the Ising energy and not by the microscopic time scale. 
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The kinetic Ising model is a very simple but effective model to study non-equilibrium 
situations. The model is based on the Ising model which describes the interaction of a set 
of spins with values {(Xj = ±1} at each lattice site i. Starting from an arbitrary initial state, 
the simplest dynamics consist of a single spin-flip process (Xj —cXi which is realized with 
a certain transition rate Wi{ai). This rate is included into a Master equation, which is an 
equation of motion for the single-time probability p{n,t). In our case the configuration n 
consists of the set of all spins in d dimensions. This is usually referred to as the Glauber 
model Q]. The problem is to find out an analytical expression for the transition rate. In 
order to ensure that the system eventually relaxes to an equilibrium state, one imposes 



the principle of detailed balance, and 
principle. One choice is proposed in 



;he transition rate is chosen in accordance with that 



Wi) = — 1 — CTi tanh — 



(1) 



Here the temperature T is given in terms of the Boltzmann constant and Ei is the local 
energy of the Ising model 

Ei = hi + '^ Jij(Tj , (2) 

where means the sum over all nearest neighbors of lattice site i. The local energy arises 
from the Hamiltonian given by 

H = hiOi Jij(Ti(Tj . (3) 

The quantity hi is an external field, and the summation in the interacting part goes over 
all pairs of nearest neighbor spins. Notice that the choice of the transition rate in Eq. 
is not unique. Moreover, the interaction energy and the coupling to the heat bath with 
temperature T is only incorporated into the Master equation via the principle of detailed 
balance. Let us stress that the different choices for the transition probability is discussed 
in the context of Monte Carlo methods j^. The goal of the present note is to include the 
energy functional and the temperature directly into the Master equation from the beginning. 
To that aim we use a mapping of the Master equation onto a dynamic equation in terms of 



seco.. ,,a...ed op.ato.4HaHa, fo. a _t .View see H. 

Let us start from a general Master equation written in the form 



dtp{n, t) = [w{n I n') p{n', t) — w{n' \ n) p{n, t) ] = L{n, n') p{n', t) . (4) 



Here p{n, t) is the probability that a certain configuration n is reahzed at time t and w{n \ n') 
plays the role of the transition probability per unit time from configuration n to n' . In our 
case the configuration n is given by the orientation of the set of spins. The principle of 
detailed balance means that the stationary distribution Ps{n) fulfills 

(5) 
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In case the static properties of the system are governed by a Hamiltonian let us make the 
following ansatz 

w{n' I n) = exp[-(3H{n')/2] V{n \ n) exp[l3H{n)/2] . (6) 

This ansatz is motivated by the conventional Arrhenius ansatz for transition rates, see also 

the result obtained in Eq. (fTT^ and the comment made there. Further let us remark that 

we have assumed the validity of the commutator relation [H{n), H{n')] = 0. In Eq. ^ the 

parameter j3 is an arbitrary one. In case the system is coupled to a heat bath we will identify 

f] with the inverse temperature T in units of the Boltzmann constant. Then the stationary 

condition Eq. can be rewritten as 

Ps{n') V{n' I n) 

= T/^ I f\ exp{f3[H{n) - H{n . (7) 
Ps{n) V{n I n') 

Inserting Eq. (jH)) in the Master equation (jH) we get 

dtp{n,t) = J2[^M-f3H{n)/2)V{n\n') exp{(3H{n')/2)p{n',t) 

n' 

- exp{-(3H{n')/2)V{n' \ n) exp{f3H{n)/2) p{n, t) ] . (8) 
Now the further aim is to rewrite Eq. JHI) using second quantized operators. Following Refs. 

innnnn h« & h ^ 

^, |5|, la |8|, |10|, lUl , the probability distribution p{n,t) can be related to a state vector \F(t)) 
in a Fock-space according to p{n,t) = {n\F(t)) and \F(t)) = J2nP(^'^) l'^)' respectively, 
where the basic vectors \n) can be expressed by second quantized operators. Using this rep- 
resentation, the underlying Master equation can be transformed into an equivalent evolution 
equation in a Fock-space, written in the form: 

d,\F{t)) = L\F{t)) . (9) 

The dynamical matrix elements L{n,n') within the Master equation are mapped onto the 
operator L = L{a,a^), where a and a"^ are the annihilation and creation operators, respec- 
tively. Here the matrix elements of the operator L(a, a^) coincide with the matrix elements 



L(n, n'). Oriffinally, this transformation had been apphed for the Bose case with unrestricted 

n n n 

occupation numbers j^|5|, |6|]. Here, we consider the case of restricted occupation numbers 
0, y, Q, Q] • In order to preserve the restriction of the occupation number in the underlying 
dynamical equations, the commutation rules of the operators a and are chosen as Pauli 
operators 0, sl. Iisj]: 



[oi, a\] = 6ij{l - 2alai) , [a^, aj] = [a], aj] = , = (a|) 



U2 



0. 



(10) 



The relation to the spin variable is = 1 — 2al at. In case of a single spin-flip process the 
evolution operator L reads 



L = ^[A(l-al)a, + 7(l-a,)al 



(11) 



where A and 7 are the temperature dependent transition rates, the determination of those is 
beyond the scope of the present approach. The flip-rates are assumed in accordance to the 
principle of detailed balance manifested in Eq. Obviously, the transition rates should 
depend on the details of the mutual interaction of the spins. Therefore the evolution operator 
should be extended by including the temperature and the interaction. In accordance with 
Eq. (jSJ we propose the following generalization 



L = [(1 - 4) exp{-pH/2)aiexp{pH/2) + (1 - a,) exp{~f3H/2)al exp{(3H/2) 



(12) 



Here, k is a parameter which flxes the time scale of the flip-process and H is the Hamiltonian 
for the underlying interaction given by Eq. (jH]). Because of the relation between ai and the 
annihilation and creation operators the evolution operator can be rewritten in terms of these 
operators as 



with Ei = hi + J(0) - 2 ^ Jija^aj ; J(0) = ^ Jij . 



(13) 



Notice that the last relation is derived only by using the algebraic properties of the operators. 
In case of vanishing mutual interaction, i.e. J = 0, one observes that the last relation is 
equivalent to the conventional Arrhenius ansatz. To that aim the comparison of Eqs. (fTTj) 
and (fT^ yields to the identification A = fi;exp( h/T) and 7 = Kexp{—h/T) in according to 



the Arrhenius ansatz. Our quantum approach yields therefore a the possibihty to formulate 
the Arrhenius ansatz in a more formal and mathematical manner. 

To proceed further we follow Doi and calculate the average of an arbitrary physical 
quantity B{n) by using the average of the corresponding operator B = \n) B{n) {n\ via 
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im) = Y.p{n,t)B{n) = {s \B\ F{t)) . (14) 

n 

Here we have used the projection state (s| = Yin (^1' which is realized only for spin 1/2 
fermions in such a simple form. The normalization condition for the probability density is 
included in the condition {s\F{t)) = 1 with the consequence that the evolution operator 
fulfills always the relation {s\L = 0. In the present case the averaged spin- variable obeys 

2!:|<--)-(-K?)>-f-^°*(?)>- 

Based upon the conventional Master equation Eq. Q the averaged spin satisfies 

d 

— {ai) = -2{aiWi{(Xi)), (16) 

where due to Glauber jll] or Suzuki et al. the single transition rate Wi{ai) = 

w{ai, . . . , (Ji , . . . , (Tat I 0"!, . . . , —(Ti , . . . , (Tn) is heuristically chosen in Eq. ((T}. This form 
of the transition rate is not uniquely determined by the principle of detailed balance. In our 
approach we can directly find the transition rate by applying the relation 

EiC.; \ , / E.; \ . . I E,; 



Oi exp ( — \ = Oi cosh ( y 1 - sinh ( — 



Using this relation and Eqs. (jlSpifij) we get immediately 



Wi{oi) = Kexp{-^^) = /tcosh 



1 — cr,- tanh — 



(17) 



T 

The transition rate is related to the heuristic one by assuming that the time scale a in Eq. 
is controlled by the spin configuration as well as the temperature. Now let us discuss the 
transition rate obtained by Eq. ()17|) in detail, especially for the case of a vanishing external 
field. In that case the local energy is reduced to Ei = Jij^j ' where the summation is 
performed over all the z nearest neighbors of the lattice site i. In the high temperature limit 
both transition probabilities wf and Wi does not distinguish. In that case the transition rate 
is obviously only determined by the microscopic time scale a or Both probabilities are 
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independent on the spin configuration. In the low temperature hmit T <^ Ei the situation 
is completely different. Let us firstly assume that the spin at site i is directed upwards, i.e. 
(7i — l. Then one has to distinguish two cases: 

(i) Ei < 0: In this realization the majority of nearest neighbor spins is not adapted to 
the preferred upward direction of the spin cxj. With other words, the local spin configuration 
around CTj is unfavorable. As a consequence the transition probabilities behave like 

lirn (cTj = 1) = fi;exp(| E^ \ /T) 

lim wf ((7i = 1) = - . (18) 

T^o a 

While in the Glauber transition rate the microscopic time scale plays the role of a lower 
cut-off, i.e. spin-flips are possible within this time scale, in our realization the transition 
rate is controlled by the by the energy E^ of the Ising model. The rate increases drastically 
and leads to an immediate flip-process of the "wrong" spins. The time scale Tj for a flip is 
of the order 

n ^ K-^ exp(- I Ei I /T) . 

The spins in the energetically unfavorable direction perform the flips with a very high rate 
in a quite short time interval. 

(ii) Ei > 0: In that case the majority of the spins around lattice site i are already adapted 
and the transition rate tends to zero according to 

Wi{a^ = l) = Kexp{-E,/T). (19) 

A similar behavior is also observed within the approximation due to Glauber. 

In the second case we consider the downward orientation ai — —1. As above we study the 

two cases of positive and negative local energy Ei. For ii^j < we get in the low temperature 

regime 

^rawi {ai = -1) = Kexp(- \ Ei \ /T) ^ , (20) 

The transition rate tends to zero, i.e. all spin flips are suppressed completely. In the opposite 

case Ei > it results 

lim (cTj = -1) = Kexp{Ei/T) , (21) 
which supports a very high flip rate, whereas the Glauber rate remains simply constant. 
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Summarizing our brief report, we have revisited the well established kinetic Ising model in 
terms of second quantized operators. The coupling to a heat bath at temperature T and 
the underlying interaction are included similar to the Heisenberg-picture of operators. As 
the result we get an exact expression for the transition rate which is likewise in accordance 
to detailed balance. The differences to the conventional Glauber model consists in the low 
temperature regime. Here the transition rate is controlled by the Ising energy and not by 
the microscopic time scale. For low temperatures the transition rate is better adapted to 
the physical situation in mind. 
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